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Summary. In the last decade it has been proven that the standard spin-wave 
theory was able to provide accurate zero-temperature results for a number of low- 
dimensional Heisenberg spin systems. In this chapter we introduce the main ingre- 
dients of the spin-wave technique using as a working model the two-leg mixed-spin 
ferrimagnetic ladder and the Dyson-Maleev boson formalism up to second order in 
the spin-wave interaction. In the remainder, we survey typical applications in low- 
space dimensionality as well as some recent modifications of the theory admitting 
a quantitative analysis in magnetically disordered phases. The presented spin-wave 
results are compared with available numerical estimates. 



1 Introduction 

The spin-wave theory is probably one of the most powerful tools ever used in 
the theory of magnetism. Originally proposed by Bloch [1, 2] and Holstein and 
Primakoff [3] as a theory of the ferromagnetic state, it was later extended for the 
antiferromagnetic Neel state by Anderson [4], Kubo [5], and Oguchi [6]. Dyson's 
profound analysis of spin-wave interactions [7, 8] demonstrated that spin waves may 
be used to obtain asymptotic expansions for the thermodynamic functions of the 
Heisenberg ferromagnet at low temperatures. Dyson's method was generalized by 
Harris et al. [9] to calculate in a systematic way spin-spin correlations, spin-wave 
damping, and various thermodynamic properties of antiferromagnetic insulators. 

It should be noticed that the basis of the spin-wave theory (SWT) for antiferro- 
magnets is much less established than for ferromagnets. The Dyson-Maleev trans- 
formation [10] gives a correspondence between any operator defined on the Hilbert 
space of the spin system and an operator on the boson Hilbert space. Evaluating 
the required averages for the Bose system, we necessarily make two approximations. 
First, we expand these quantities, by using a perturbation formalism in which the 
unperturbed Hamiltonian is quadratic in boson operators and the perturbation is 
the remaining quartic interaction. Second, we neglect the projection operator in the 
averages, which takes into account the so-called kinematic interactions by canceling 
the boson states with more than 2S bosons per lattice site, S being the spin quan- 
tum number of the lattice spin. In the ferromagnetic case, Dyson has argued that 
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these approximations would lead to results which are asymptotically correct at low 
temperatures (T) to all orders in T. In the antiferromagnotic case, the situation is 
less settled due to the zero-point motion, i.e. quantum spin fluctuations in the Neel 
state. In principle, one may suspect that there are errors in the perturbation theory 
even at zero T. The same problem appears in the Holstein-Primakoff formalism 
[3]. We refer the interested reader to the original papers cited above as well as to 
the monographs [11, 12, 13] for details concerning this problem. In principle, the 
spin-wave approach is less efi^ective for low-dimensional quantum spin systems, as 
quantum spin fluctuations typically increase in reduced space dimensions (D) and 
for small spin quantum numbers S. Moreover, since at finite T thermal fluctuations 
completely destroy the magnetic long-range order in ID and 2D Heisenberg models 
with isotropic short-range interactions [14], in such cases the conventional SWT 
completely fails. 

In view of the mentioned drawbacks of SWT, it seems surprising that for the 

last decade the standard spin-wave approach has been found to give very accu- 
rate description of the zero-temperature physics of a number of low- dimensional 
spin models, the best example being the S = ^ Heisenberg antiferromagnet on a 
square lattice [15]. Probably, another good example is the mixed-spin Heisenberg 
chain describing a large class of recently synthesized quasi-lD molecular magnets 
[16] (cf. Chap. 4). The following analysis reveals some common features of these 
examples, the most important being the weakness (in a sense) of spin-wave inter- 
actions. Fortunately, in low-space dimensions many numerical techniques - such 
as the quantum Monte Carlo method (QMC), the exact numerical diagonalization 
(ED), and the density-matrix renormalization group method (DMRG) - are more 
effective, so that the discussed drawbacks of the spin-wave analysis may be partially 
reduced by a direct combination with numerical methods. 

A goal of the present review is to summarize typical applications and some re- 
cent developments of the spin-wave approach related to low-dimensional quantum 
spin systems. The spin-wave technique is presented in the following Sect., using the 
mixed-spin Heisenberg ladder as a working model and the Dyson-Maleev boson 
formalism. Due to the asymptotic character of spin-wave series, the calculation up 
to second order in the spin-wave interaction is a reasonable approximation for most 
of the applications at zero T. As fax as at this level perturbative corrections can 
easily be calculated in the framework of the Rayleigli Sclirodinger theory, we will 
not consider in detail perturbation techniques based on niagnon Green's functions 
[9, 17]. Typical applications of the spin- wave formalism in low-dimensional spin 
systems are presented in Sects. 3 and 4. In particular. Sect. 3 involves an analy- 
sis of the parameters of the quantum ferrimagnetic phase in mixed-spin quasi- ID 
models, such as the (,si,,S2) Heisenberg chain. The SWT results are compared with 
available DMRG and ED numerical estimates. Section 4 collects basic SWT re- 
sults concerning 2D Heisenberg antiferromagnets. Some recent modifications of the 
SWT - admitting a quantitative analysis in magnetically disordered phases - are 
presented in Sect. 5. Section 6 contains concluding remarks. 

2 Dyson— Maleev Formalism 

In this Sect, we describe the formal apparatus of the SWT. We choose as a working 
model the mixed-spin Heisenberg ladder (Fig. 1) defined by the Hamiltonian 
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JV JV 
'H = ' ^"+1 + ■ «n+ll + Jl. Sn- CTn, (1) 

n=l n=l 

where the index n (= 1, • • • , A'^) labels the rungs of the ladder, and N is an even 

integer. The ladder is composed of two typos of spins (s„,cr„) characterized by 
the spin quantum numbers si and S2 (si > S2): s„^ = Ti^si{si + 1) and ffn^ = 
Tt? S2{s2 + 1). In the following text we use the notation Va = si/s2 > 1, and set 
Ti = 1 and oo = 1, oo being the lattice spacing along the ladder. 



n+1 



Fig. 1. Mixed-spin Heisenberg ladder composed of two types of site spins. The 
arrows show one of the classical ground states for J_l > 0, defined by the orientation 
of the ferromagnetic moment M. = X^„^i(sn-|-<Tn)- The intrachain coupling J = 1. 



It is worth noticing that the model (1) is not purely academic. For instance, 
recently published experimental work on bimetallic quasi-lD molecular magnets 
(cf. Chap. 4) implies that the magnetic properties of these mixed-spin compounds 
are basically described by the Heisenberg spin model with antiferromagnetically 
coupled nearest-neighbor localized spins. The ladder structure in Fig. 1 reproduces, 
in particular, arrangements of the Mn {s\ = |) and the Cu (s2 = \) magnetic 
atoms along the a axis in the compounds MnCu(pbaOH)(H20)3 (pbaOH — 2~ 
hydroxy-1 ,3-propylenebisoxamato) [18] . 



2.1 Clcissical Reference State 

The first step in constructing a spin-wave expansion is to find the lowest-energy clas- 
sical spin configurations of the related classical model. As a rule, this is a straight- 
forward task, apart from some magnetic models with competing interactions which 
may exhibit complicated non-collinear spin states (see, e.g. [19]). Another serious 
problem at this stage may be related to a macroscopic degeneracy of the classical 
ground state, a typical example being the Heisenberg model on a kagome lattice (cf. 
Chap. 3) which exhibits a magnetically disordered ground state. Further analysis 
of the problem involves quantum fluctuations and the so-called order-from-disorder 
phenomenon [20, 21]. 

Turning to our model (1), it is easy to see that the required reference state 
for Jx > is a ferrimagnetic spin configuration where the s„ spins are oriented 
in a given direction, and the a-n spins point in the opposite direction (see Fig. 1). 
The state is degenerate under arbitrary rotations (as a whole) in the spin space. 
One may pick up a reference state by introducing a small staggered field, say, for 
the Sn spins. We can actually get more information even in the quantum case, by 
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using the Lieb-Mattis theorem for bipartite lattices [22]. First, the theorem predicts 

that the quantum ground state belongs to a subspace with the total-spin quantum 
number [Si — S2)N, i.e. for J± > the system has a ferrimagnetic ground state 
characterized by the ferromagnetic moment per site Mo = (si — S2)/2. Second, the 
theorem states that the energies of the ground states E{St) characterized by the 
total-spin quantum numbers St > N{si — S2) are arranged as follows 

E{St + 1) > E{St) ■ (2) 

Notice that the classical and quantum ferrimagnetic ground states have one and 
the same magnetization Mo, but otherwise they are different because the clas- 
sical ground state is not an eigenstate of the quantum model (1). The quan- 
tum ferrimagnetic state is [2A'^(si — S2) -I- l]-fold degenerate, since the z com- 
ponent of the total spin - being a good quantum number - takes the values 
—N{si — S2), —N{si — S2) + 1, • • • , N{si — 82). This quantum magnetic phase may 
also be characterized by the following sublattice magnetizations 

^ JV 1 ^ 

mA = -^^(s„) ms = -^^(<Tn), (3) 

n=l n=l 

where the symbol (• • •) means a quantum-mechanical average over the ground state. 
We shall later see that quantum spin fluctuations reduce the classical sublattice 
magnetizations si and S2, but the magnetic long-range order is preserved, i.e. 
TTlA, niB / 0. 

In the region J± < the situation is different, i.e. the lowest-energy spin con- 
figuration is the Neel antiferromagnetic state based on the composite rung spins 
si + S2- Now the Lieb Mattis theorem predicts that the exact quantum ground 
state is a spin-singlet state, i.e. St ~ and Mo — 0. Therefore, it may be gener- 
ally expected a magnetically disordered phase, t.e. mA,mB = 0, as the isotropic 
Heisenberg model (1) is defined on a bipartite ID lattice (see, e.g. [23]). In terms 
of the SWT this would mean that the classical antiferromagnetic state is swept out 
by quantum fluctuations, so that the concept of the spin-wave expansion does not 
work at all. 



2.2 Boson Hamiltonian 

Now we describe the second step in constructing the spin-wave expansion, t.e. the 
transformation of (1) to a boson Hamiltonian. The most popular boson represen- 
tation of spin operators has been suggested by Holstein and Primakoff [3]. Other 
useful representations have been devised by Schwinger [24], Maleev [10], Villain 
[25], and Goldhirsch [26, 27]. 

We start by defining the Holstein-PrimakofF representation for the spins s„ 
(n=l,...,7V): 



s„ = V2si Y 1 ^j— a„ , s„ = V2si a^^W 1 ^j— , s„ = si-a>„ , (4) 

where = ± and si is the spin quantum number. a„ and are annihilation 
and creation boson operators satisfying the commutation relations 
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[a„, a]ji] = 5„m, [a-n, am] = [at, oL] = . (5) 

Using the last equations, it is easy to show that the operators defined by (4) satisfy 
the commutation relations for spin operators 

and the equation = si (si + 1). The operators a„ and at act in the infinite- 
dimensional boson Hilbert space spanned by the orthonormal bcisis states 

|ni,n2,...,niv)= ^ / f , —Ji' |0) , (7) 

where rii (= 0, 1, . . . , oo) is the occupation number of site i. The reference vacuum 
state |0) is defined by the relations ai\0) = (for V i). 

It is possible to rationalize the square roots in (4) by the Maleev similarity 
transformation 

t \ 1/2 / t \ -1/2 

This transformation is not unitary, but preserves the number operator atari as 
well as the commutation relations (5) within the physically relevant Hilbert space 
{rii < 2si for V i). Applying the last transformation to (4), we get the Dyson-Maleev 
boson representation 

= V2si{l- ata„/2si)an , s~=V2siat, = si - atan ■ (9) 

Note that the operators s^ in this representation are not Hcrmitian conjugate in the 
boson space (7) so that in the general case they will generate non-Hermitian Hamil- 
tonians. Treatment of such Hamiltonians requires some care, but it seems that - at 
least up to second order in the spin-wave interaction - this does not cause serious 
problems. More problematic is the relation between physical and unphysical states. 
The latter appear in the exact Holstein-Primakoff representation as well, as any 
actual calculation requires truncation of the asymptotic square-root series. Dyson's 
method [7] eliminates the unphysical boson states by a projection operator giving 
zero on these states. In practice, however, we are enforced to eliminate this operar 
tor. As already mentioned, this is the basic approximation of SWT. As a whole, the 
Dyson Maleev formalism has many advantages if one needs to go beyond the linear 
spin-wave theory (LSWT) within a perturbation scheme. This is because the inter- 
actions between spin waves are better handled so that the unphysical singularities 
caused by the long- wavelength spin waves cancel out. 

To continue, we write a representation similar to (9) for the spins (t„, by using 
a new set of boson fields (6„, n = 1, . . . , TV): 

a:^ = V2s^bt{l-btbn/'2s2), ff^=V2s^bn, cr^ = -S2 + btbn ■ (10) 

b„ and bt satisfy the same commutation relations (5), and are supposed to commute 
with the set of a bosons. Here the reference state is chosen in the opposite direction, 
in accord with the classical spin configuration in Fig. 1. 

Using (9) and (10), we can find the boson image of any function of spin op- 
erators. In particular, we are interested in the boson representation of the spin 
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Hamiltonian (1), which we denote by Hb- For the purposes of SWT, it is instruc- 
tive to express 7is in terms of the Fourier transforms ak and hk of the boson 
operators a-n and &„, by using the unitary Fourier transformations 

fc k 



and the identity 



1 i(k-k')n s 



where 



=1 

It may be verified that this transformation is canonical, by showing that the new 
operators ak and bk obey a set of commutation relations identical to (5). The wave 
vectors k in the last expressions are defined in the first Brillouin zone: 

fc- <- 2 2 +^'---' 2 ' • 

Notice that the rung spins {sn,a„) in Fig. 1 compose the n-th magnetic (and 
lattice) elementary cell: this may be easily observed by interchanging the site spins 
of every (say) even rung in Fig. 1. 

We leave the Fourier transformation oiTi.B as an exercise, and directly present 
the result in terms of the new operators Ofc and bk'- 

■Kb = -2'rorsS^ +no + Vdm , (12) 
7^0 = [70 {alak + r«&l6fc) + Vrljk {albl + akbk)] , (13) 

k 

Vdm = --^^ 512(271-44026! 64 + V»^7i+2-4a3b2&I 64 + -^744020164 j • (14) 
1-4 ^ V s / 

Here jk = J± /2 + cos k (70 ~ J±/2 + 1), = A{ki + fe2 — fca — fc4) is the Kronecker 
function, and we have introduced the abbreviations (fci, fe2, fea, fe4) = (1,2,3,4) for 
the wave vectors. 

In a standard spin- wave expansion, 1/si and l/s2 are treated as small parame- 
ters, whereas the parameter r-g may be considered as a fixed number of order unity. 
In such a perturbation scheme, it is convenient to set 1/5 = I/S2 and use 1/5 as 
a small parameter. Thus, the first term in (12) - the classical ground-state energy 
- is proportional to S^, the LSWT Hamiltonian Tio is multiplied by S, and the 
spin-wave interaction term Vdm ha^s the order 0(1). We shall follow a perturbation 
scheme where the diagonal terms of V^mj i-e- terms proportional to the occupation- 
number operators a^ak and blbk, are treated together with Ho as a zeroth-order 
Hamiltonian, whereas the rest of Vom is taken as a perturbation [9]. This is a more 
generic approach because for some reasons the spin-wave interactions may be weak 
even in the extreme quantum systems with 1/S = 2. 
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2.3 Quasiparticle Representation 

In the next step, we diagonalize the quadratic Hamiltonian Ho, by using the Bo- 
goliubov canonical transformation to quasiparticle boson operators (a^ and f3k) 
[3]: 

ak = Uk{ak - XkPl) , hk = UkiPk - Xka\) , ul{l-x'i) = l. (15) 

It is a simple exercise to find the transformation parameters Uk and Xk from 
the condition which eliminates the off-diagonal terms akl3k appearing in Ho after 
the transformation (15). The result reads 



where 

In some applications, the quadratic Hamiltonian Ho may include additional ferro- 
magnetic bilinear terms (such as a\.hk) so that the actual diagonalization is more 
involved due to the increased number of parameters (16). Some diagonalization tech- 
niques for systems with large number of boson operators are presented in [11, 28]. 

A quasiparticle representation of the quartic terms (14) requires more techni- 
cal work. As mentioned above, it is instructive to pick up the quadratic diagonal 
terms in Vj^m and to treat them together with Ho as a zeroth-order approxima- 
tion. A simple way to do this is based on the presentation of Vjjm ^ a sum of 
normal-ordered products of boson quasiparticle operators. Apart from a constant, 
the resulting expression for Vj^m contains diagonal and off-diagonal quadratic op- 
erator terms, and normal-ordered quartic operator terms. We leave as an exercise 
this simple but somewhat cumbersome procedure and give the final result for Hb 
expressed in terms of the quasiparticle boson operators and 

Hb=Eo+Hd + \V , V = V2 + Vdm, Asl. (18) 

Here Eo is the ground-state energy of the ferrimagnetic state calculated up to 
the order 0(1) in the standard 1/5 expansion: 

^ = -27or.S^ - 7o(l + r.) ||l - 1 ^ j 5 + ei + O (i) , (19) 
where ei = — 2(c? -|- c|) — J±(c? -|- cl) — (2c2 -|- J±C3)ci(rs -|- l)r'^^^^ and 

k k k 

Hd is the quadratic Hamiltonian resulting from Ho and the diagonal terms 
picked up from (14): 



Hn=2SY,[4"^»Wk+Jk'^0lPk 



(21) 



where up to O (1/5) the dressed dispersions read 
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„<..«.,(i±i..,r^)4 + c(J,) (2.) 

where 5^ = {gkVk - do)£'^^^'^ /2 ± (r^ - l)(2c2 + C3 Jx)r7^^^/2, gfe = 2ci(rs + 
1)7^''^^''^ + 4c2 cosfc + 2c3Jx, do = 4ci7o + (r^ + l)(2c2 + J±C3)rs^^^. 

The functions Wj."''^^ without 0{1/S) corrections will be referred to as bare 
dispersions. 

Finally, the quasiparticle interaction V includes two different terms, i.e. the 
two-boson interaction 

V^ = Y1 [Vk'^ll^l + Vk-akl3k] (23) 

k 

defined by the vertex functions 

and the quartic Dyson-Malccv interaction 

VdM = '^12 |^V/2.34Q!|a2a3a4 + 2V'/2^.34a|/?2a3a4 + 2V/2%4Q!ia2/3304 

1-4 

+ 4^^340+03/31/32 + 2F/2%/3la3/32/3i + 2F/2%/3l/3la+/3i 

(25) 

defined by the vertex functions V^2.34, i = 1, . . . , 9. We have adopted the symmetric 
form of vertex functions used in [17]. The explicit form of ^^2.34 depends on the 
concrete model. For the ladder model (1), the vertex functions may be obtained 
from those of the Heisenberg ferrimagnetic chain [29] , using the formal substitution 
cos k I — > cos k + J±/2. 

In the following we shall treat the spin-wave interaction V as a small pertur- 
bation to the diagonal Hamiltonian Eq + Hd- To restore the stamdaxd 1/S series, 
one should (i) use bare dispersion functions, and (ii) resume the series in powers of 
1/5. 



3 Spin Wave Analysis of Quasi-ID Ferrimagnets 

In this Sect, we analyze the magnon spectrum and basic parameters of the quantum 
ferrimagnetic phase of the model (1), by using the developed spin-wave formalism 
and the Rayleigh-Schrodinger perturbation theory up to second order in A. The 
SWT results are compared with available DMRG and ED numerical estimates. 

3.1 Linear Spin Wave Approximations 

In a standard linear spin-wave approximation we consider only the first two terms 
in (12), and discard Vqm as a next-order term iu 1/5. This corresponds to the 
first two terms in the expression for the ground-state energy (19), and to the first 
term in the expression for the quasiparticle dispersions (22). As a matter of fact, 
by using the normal-ordering procedure, we have already got even the next-order 
terms of the expansions in 1/5 for these quantities. 
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Magnon Excitation Spectrum 

The quadratic Hamiltonian Hd defines two branches of spin- wave excitations (a 
and /? magnons) described by the dispersion functions cj^" "^^ in the first Brillouin 
zone — TT < fe < TT (see Fig. 2). The excited states a]^|0) belong to the 

subspace characterized by the quantum number S J = St — 1 {S^ = St + 1), where 
St = (si — S2)N . In the long wavelength limit A; <C 1, the energies of a magnons 
E^^^ have the Landau-Lifshitz form 

E^^^^2SJ-^ = j^^k'+0{k'), (26) 

where Qs is the spin stiffness constant [30]. This form of the Goldstone modes is 
typical for Heisenberg ferromagnets, and reflects the fact that the order parameter, 
i.e. the ferromagnetic moment, is itself a constant of the motion. 




Fig. 2. Magnon excitation spectrum of the mixed-spin ladder (,si,S2) = (1, \) for 
interchain couplings J± = 0.1 and J± = 1. The dashed lines display the energy of 
P magnons e\^^ related to the Hamiltonian T-io- The solid lines show the magnon 
spectra as obtained from the second-order approximation in V . The energy of a 
magnons related to (22) is not displayed, as it closely follows the respective solid 
lines. The symbols indicate ED numerical results. The Figure is taken from [31]. 



The spin stiffness constant as well as Mo play a basic role in the low- 
temperature thermodynamics [32]. The parameter may be obtained from the 
Landau-Lifshitz relation and (22): 

" =^-lic.'-^^^^^0(l,\. (27) 



exhibits 

in the vicinity of vr it reads 



The function -E^"' exhibits an additional minimum at the zone boundary, so that 



^(a) ^ ^(a) ^ ^^^^^ _ ^^2 ^ ^-28) 
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Here is the excitation gap at the zone boundary. In the Umit J± — > 0, the 

excitation gap zli-"' (cx J±) goes to zero. For ferromagnetic couplings J± < 0, the 
k — n mode becomes unstable and produces global instability of the ferrimagnctic 
phase. 

The function E^f'' = 2Su)\f'' may be characterized by the spectral gaps A^^^ 
(at fc = 0) and (at fc = tt). The expression for A[P reads 

For the (si,S2) — (1,|) chain (J^ = 0), the last equations give the results 
Qa/2siS2 = 0.761 and A^^^ = 1.676, to be compared with the results Qs/2siS2 = 1 
and A\!^^ = 1 obtained in a standard linear approximation using the Hamiltonian 
Tio [33, 34]. A comparison with the numerical QMC result A\^^ = 1.759 [35] clearly 
demonstrates the importance of the 1/5 corrections to the dispersion functions (22) 
in the extreme quantum limit. 

Summarizing, it may be stated that the linear approximation - based on the 
quadratic Hamiltonian Tin - gives a good qualitative description of the magnon 
excitation spectrum of the model (1). The same conclusion is valid for the ground- 
state energy: The expression (19) has been found to produce an excellent fit to the 
numerical ED results in a large interval up to Jx = 10 [31]. 



Sublattice Magnetizations 

The on-site magnetizations tua ~ (s^) and tub = —{(^n) a^re parameters of the 
quantum ferrimagnetic phase which keep information for the long-range spin cor- 
relations. The simple LSWT results tua = si — ci and ms = S2 — ci show that 
quantum spin fluctuations reduce the classical on-site magnetizations already at 
the level of non-interacting spin waves. Tio produces the same results. The ratio 

S2 — TOS _ Cl 

may be used as a measure of the zero-point motion in the quantum ground state. 
Thus, there appears to be a well-defined scmiclassical limit S ^ oo where Ho is 
a sufficiently accurate approximation, provided 01/5 <C 1. In this connection, it 
seems surprising that the spin-wave series for the 5=1 square-lattice Heisenberg 
antiferromagnct produces the excellent result tua ~ 0.3069(2) [36] - the recent 
stochastic-series QMC estimate is 0.3070(3) [37] - in spite of the fact that in this 
case the parameter Ci/S ~ 0.393 is not small. Even more illuminating is the (1, |) 
ferrimagnetic chain: In spite of the large parameter Ci/S Ri 0.610, the second-order 
SWT gives the precise result tua = 0.79388 [38] (DMRG estimate is tua = 0.79248 
[39]). It is difficult to explain the accuracy of SWT in terms of the standard 1/S 
series. However, as will be shown below, the quasiparticle interaction V produces 
numerically small corrections to the principal zeroth-order approximation. 

In the mixed-spin model (1) there appears an important first-order correction 
to the sublattice magnetizations which is connected to the quadratic interaction 
V2. Let us go beyond the linear approximation and calculate the 0{X) correction 
to ruA- The on-site magnetization ms may be obtained from the exact relation 



(30) 
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Fig. 3. On-site magnetization (sublattice A) of the (1, i) ladder as a function of 

the interchain coupling J±. The dashed and dashed-dotted lines display the series 
results up to first order in 1/5 (bare dispersions) and V (dressed dispersions). The 
solid line shows the series result up to second order in V . The Lariczos ED results 
for ladders with A'^ = 12 rungs are denoted by open circles. The Figure is taken 
from [31]. 



ruA = si — S2 + resulting from the conservation law for the ferromagnetic 
moment. The expression of ruA in terms of quasiparticle operators reads 



(31) 



mA = s^-c,-^^[l-{alak+ - J (at + /jt/jt ) 

k 

Now we make use of the standard perturbation formula 

^Qyi) _y mv\n}{n\d\0} ^ y (n\V\0)(0\d\n) ^g^) 



giving the first-order correction in V of (O). Here O is an arbitrary operator and 
(• • •) means a quantum-mechanical average over the exact ground state. The formula 
is also valid in the case of non-Hermitian perturbations V. In our case, O is a 
quadratic operator, so that the sum in (32) is restricted to the two-boson eigenstates 
\nk) = CKfc/JfciO) of TId, k being a wave vector from the first Brillouin zone. The 
energies of these states are Ek-Eo = 2S{oji"^ + '). Finally, using the matrix 
elements 

(0\V2\nk) = , (nfc|V2|0) = V^+^ , (33) 
we get the following result for niA calculated up to first order in V: 



ruA = si — ci 



1 ^VkV^^^+vr 



4SN _ 



To find the standard 1/5 correction to tua, we have to use in (34) the bare dispersion 
functions. 
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Figure 3 shows the results for ruA, as obtained from (34) by using the bare and 
dressed dispersion functions (22). It is seen that the expansion in 1/S gives a small 
(unexpected) decrease of niA in the vicinity of J± — 0, whereas the expansion in 
V produces a correct qualitative result in this limit. The indicated problem of the 
standard 1/S series probably results from enhanced fluctuations of the individual 
chain magnetizations about the common quantization axis. Indeed, at the special 
point Jx = the classical ground state acquires an additional degeneracy under 
independent rotations of the chain ferromagnetic moments. Thus, the quartic di- 
agonal interaction - included in Wo - seems to stabilize the common quantization 
axis connected to the global ferromagnetic moment. We have an example where the 
expansion in powers of V gives better results. 



Antiferromagnetic Chain 



It is instructive to consider the antiferromagnetic chain as a special case (si = 
S2, J± = 0) of the mixed-spin model (1). After some algebra, from (19) and (22) we 
find the following simplified expressions for the ground-state energy (per site) 



eo 



1 + 



and the magnon spectrum 



,(«,/3) 



2S 



1 

2S 



1 - 



Isinfcj +0 



(35) 



(36) 



of the antiferromagnetic chain. For S = ^, the standard LSWT gives the result 
eo = -0.4317 which is close to Hulthen's exact resuh - In 2 -I- 1/4 w -0.443147 
[40]. It is an illuminating agreement, as the theory might have been expected to 
fail for magnetically disordered states. Notice, however, that the next-order ap- 
proximation, i.e. Co = —0.4647, does not improve the SWT result. This indicates 
a poor convergence of the 1/S expansion. We can also check the series for 5 = |, 
by using the numerical result eo = —2.82833(1) [41] based on DMRG estimates 
for finite systems and the finite-size corrections to the energy, as derived from the 
Wess-Zumino-Witten theory [42]. The first two terms in the series (35) for 5 = | 
give the result eo = —2.79507. In this case, an inclusion of the next-order term in 

(35) produces the precise SWT result eo = -2.82808. Thus, aheady for 5* = | the 
spin-wave series shows a good convergence. 

Turning to the magnon spectrum (36) , we find that for S = ^ SWT qualitatively 
reproduces Des Cloizeaux and Pearson's exact result for the one-magnon triplet 
excitation spectrum Ek = ^\smk\ [43]. It is interesting that the 1/S correction 
in (36) improves the standard LSWT result for the spin-wave velocity (c = 1) to 
the value c = 1.3634: the exact result is c = 7r/2 w 1.5708. The magnon spectrum 

(36) is doubly degenerate and has the rclativistic form Ek ~ c\k\ {c\tt — k\) near 
the point fc = (fe = tt), to be compared with the rigorous result where the spin- 
wave states, being eigenstates of spin 1, axe triply degenerate. Long-wavelength spin 
waves correspond to states where all regions are locally in a Neel ground state but 
the direction of the sublattice magnetization makes long-wavelength rotations. 

Using (20) and (30), we find the following expression for the on-site magnetiza- 
tion in the antiferromagnetic chain 
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- = ^-'^^-^+^-2fE^ = — ■ (37) 

Wo see that in ID the quantum correction is divergent at small wave vectors already 
in the leading LSWT approximation, no matter how large is 5. This indicates that 
the Neel state is destabilized by quantum fluctuations, so that the concept of spin- 
wave expansion fails. 

Finally, it is instructive to calculate the long-wavelength behavior of the correla- 
tion function (sn-CTn+x) ■ Using the Dyson Maleev representation and (15), one finds 

{s„-(7„+j:) = -S'^+2S{a„hn+x)-\ where (a„6„+a,) = -(1/2/^) ^^(cos fe/|smA:|) exp(ifex). 

Thus, in the limit 1 one obtains 



(38) 



This indicates that in the semiclassical limit S ^ oo the antiferromagnetic chain 
is ordered at exponentially large scales ^ ~ aoexp(7rS') [44]. Here we have restored 
the lattice spacing oo. 



3.2 Spin Wave Interactions 

We have already discussed some effects of the quasiparticle interaction V , by cal- 
culating the first-order correction to the sublattice magnetizations twa and tub- 
Notice that 0(X) corrections to the ground-state energy (19) as well as to the dis- 
persion functions (22) do not appear. Indeed, it is easy to see that the corresponding 
matrix elements (0|T/|0) and {nk\V\nk) {\nk) = al|0), or vanish as a result 

of the normal ordering of V . It will be shown below that the 0{}?) corrections 
lead to further improvement of the spin-wave results. To that end, we consider two 
examples, i.e. the ground-state energy Ea and the dispersion function ii^^"'. The 
reader is referred to the original literature for similar calculations concerning the 
parameters vriA, Qs [31], and A^^^ [29]. 

The calculations may be performed within the standard perturbation formula 

^(2) ^ y {i\V\j){j\V\i) (39) 

giving the second-order correction in V to the eigenvalue Ei of the eigenstate \i) of a 
non-perturbed Hamiltonian. In our case, the zeroth-order Hamiltonian is Eq + Hd, 
and the perturbation V is given by (12). The sum in (39) runs over the eigenstates 
ofTYc. 



Second-Order Corrections to Eq 

We consider corrections to the vacuum state \i) = |0), so that the energy Ei = Eq is 
given by (19). There are two types of O(A^) corrections to Eo which are connected 
with the interactions V2 and Vdm- 

First, we proceed with the tjuadratic interaction V2. It is easily seen that only 

the states \j) = \nk} = aj./3|.|0) produce non-zero matrix elements in (39). The 
dominator for these two-boson states reads Eq — Ek = —2S{(jj\^^ -\- w^'^'), where 
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w^"''^^ are defined by (22). Using the above results and (33), we get the following 
correction to the ground-state energy (19) coming from Vi: 

Next, we consider the Dyson-Maleev interaction Vdm- Looking at the explicit 
expression of Vdm (25), we find that only the term with the vertex function V^^^i 
(^'2-34) does not annihilate the vacuum state |0) ((0|). Thus, the sum in (39) runs 
over the four-boson eigenstates |1234) = {2\2[)-^^^al^al^/3l^plj0). The related 
matrix elements read 

(1234|l/oM|0) = -^V^^Us'tt , {0|1/z,m|1234) = -^V^ig^J?! 

Using these expressions, we find the following correction to the ground-state energy 
resulting from Vdm'- 

AS) 1/(7) 



p(2)" ^ ^ "S^A^^ ^^43;12'^12;34 

° ~ 25 ,(«),,,(<»),, ,(/3) ,(/3) ■ ^ ' 



1_4 '^'1 I "-2 I "-3 I "-4 

Notice that the second-order correction to Eo in powers of l/S is the sum of E^^ 

// 

(2) 

and E^ ' but calculated with the bare dispersion functions. 

(a) 

Second-Order Corrections to ujf. 

Now we are interested in perturbations to the one-magnon states |i) = = q|.|0). 
The calculations may be performed by following the method already used for Eq. 
Since we are treating an excited eigenstate, there appear new types of corrections 
connected to the vertex functions V/2-34 '^I'd 34- These terms may be predicted, 
e.g. by drawing the diagrams shown in Fig. 4. Notice that the graphical represen- 
tation of the vertex functions in Fig. 4 is connected to the quasiparticle operator 
forms of V2 (23) and Vdm (25). The interested reader is referred to the original lit- 
erature (see, e.g. [9, 17, 45]) where this diagram technique is explained in detail. We 
leave these simple calculations as an exercise, and directly present the expression 
for the second-order corrections to Wj,™^ 



s (a) 



^ r34 / '^43;2fe '^fc2;34 4,3;2k ^ k2;34 



jV2 '"^ I ,,(«),,,(«),, ,(/3) ,,,(« ' , ,(«) , , ,(/3) , , ,(«) ,,,(«) 



2-4 



(42) 



It is interesting to note that the vertex functions V^ \ ^kp-pk^ ^/cp pfc' ^4 



(-) 1/(2) 1/(3) 1/(8) 

43;2fc' 

and 1^43.2(1 vanish at the zone center fc = ^, so that the gapless structure of cj^"' 
is preserved separately by each of the second-order corrections in (42). Thus, we 
have an example demonstrating some of the good features of the Dyson-Maleev 
formalism. 



^ Analytical properties of the vertex functions have been studied in [46] 
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v(+) v(-) 




Fig. 4. Second-order self-energy diagrams giving the corrections to the dispersion 
function ojI.°'\ Solid and dashed lines represent, respectively, the bare propagators 
for a and (5 magnons. The Figure is taken from [29]. 



3.3 Comparison with Numerical Results 

We have already presented in Figs. 2 and 3 second-order SWT results for the 
dispersion functions Wj,"'''^ and the on-site magnetization tua of the (1, i) ladder. 
The comparison shows that the SWT dispersion functions closely follow the ED 
data in the whole Brillouin zone. For instance, the SWT result for the gap Zig'^' at 
J± —0.1 differs by less than 0.5% from the ED estimate. Turning to tua, we find a 
precision higher than 0.3% in the whole interval < Jx < 3. These arc illuminating 
results, as in the considered system the perturbation parameter l/S = 2 is large. 
To understand these results, let us consider, e.g. the A series for the spectral gap 
A^ff^ of the chain [29]: 

— — 2 ^ = 1.6756A° -I- 0.1095A^ - 0.0107A^ + 0{X*) . 

2(Sl — 32) 

Although 1/5 = 2, we see that the quasiparticle interaction V introduces numeri- 
cally small corrections to the zeroth-order approximation Tin- 



Table 1. Spin-wave results for the parameters eo — Eo/N, niA, and Z\o = 
A\P /2{s-i — S2) of different (si,S2) Heisenberg chains calculated, respectively, up 
to the orders 1/5, 1/5^^, and 1/5^. The SWT results are compared with available 
DMRG estimates which are, respectively, denoted by eo, rfiA [39], and Aq [47]. 

(si,S2) eo eo tua ffiA Ao Aq 

fl,i] -1.45432 -1.45408 0.79388 0.79248 1.7744 1.76 

''1,1) -3.86321 -3.86192 1.14617 1.14427 1.6381 1.63 

f,i) -1.96699 -1.96727 1.35666 1.35742 1.4217 1.42 

'2, i) -2.47414 1.88984 1.2938 1.29 
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Finally, in Table 1 we have collected SWT results for different ferrimagnetic 
chains. It is interesting to note that oven in the extreme quantum cases (1, |) and 
(|,1), deviations from the DMRG estimates are less than 0.03% for the energy 
and 0.2% for the on-site magnetization. Moreover, it is seen that the increase of 
rs = si/s2 - keeping S2 = 5 fixed - leads to a rapid improvement of the 1/S series. 
The above results suggest that the Heisenberg ferrimagnetic chains and ladders are 
examples of low-dimensional quantum spin systems where the spin-wave approach 
is an effective theoretical tool. 



4 Applications to 2D Heisenberg Antiferromagnets 

In this Sect, wc survey recent applications of the spin- wave approach to 2D Heisen- 
berg spin systems, the emphasis being on the ground-state parameters of the square- 
and triangular-lattice Heisenberg antiferromagnets. We shall skip most of the tech- 
nical details, as the discussed spin-wave formalism actually does not depend on 
the space dimension. As already mentioned, for the last decade SWT has been 
found to produce surprisingly accurate results for the ground-state parameters of 
the square-lattice Heisenberg antiferromagnet even in the extreme quantum limit 
yS = i . Below we collect these results and compare them with recent QMC numeri- 
cal estimates. As to the triangular antiferromagnet, it seems to be a rare example of 
magnetically frustrated spin system where the spin-wave expansion is effective. In 
this case, we also give some technical details concerning the spin-wave expansion, 
as it includes some new issues resulting from the coplanar arrangement of classical 
spins. 

4.1 Square-Lattice Antiferromagnet 

The scjuarc-latticc S = ^ Heisenberg antiferromagnet being a simple and rather 
general model to describe the undoped copper-oxide materials - has received a great 
deal of interest for the last decade. Now it is widely accepted that the ground state 
of the model is characterized by antiferromagnetic long-rage order. Thus, the role 
of quantum spin fluctuations is restricted to reduction of the sublattice magnetiza- 
tion from its classical value ^ by about 39%. ^ In a seminal work by Cliakravarty, 
Halperin, and Nelson [48] - using the renormalization-group approach to study the 
quantum non- linear a model in 2 H- 1 space-time dimension - it has been shown that 
in the so-called renormalized classical regime ksT <^ ps the thermodynamic prop- 
erties of the 2D quantum Heisenberg antiferromagnet are dominated by magnon 
excitations, so that the leading and next-to-leading corrections in fcsT/ps are fully 
controlled by three physical parameters, i.e the spin stiffness constant ps,^ the spin- 
wave velocity c, and the on-site magnetization m, calculated at T = (see also [49]). 

^ Compare with the reduction of about 42% of the classical on-site magnetization 

i in the (1, i) ferrimagnetic chain (see Table 1). 

This quantity, measuring the response of the system to an infinitesimal twist of 
the spins around an axis perpendicular to the direction of the broken symmetry, 
should not be confused with the spin stiffness constant of the ferromagnetic state 
Qs connected to the Landau-Lifshitz relation (26). 
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Moreover, it has been argued that the discussed universal thermodynamic proper- 
tics appear for arbitrary ksT/ps, provided that < Ps <C J and fcsT <^ J, J being 
the nearest-neighbor exchange constant [50]. 

The quantities ps , and c appear as input parameters in the quantum non-linear 
cr model defined by the Lagrangian density 




where the vector staggered field n = n{t, x, y) satisfies the non-linear constraint 
= 1. This model may be introduced using arguments based on general grounds: 
As long as the continuous 0(3) symmetry is spontaneously broken, the symmetry 
of the problem requires that the interaction of the Goldstone modes, i.e. spin waves, 
of the system in the long-wavelength limit be described by this model regardless 
of the details of the macroscopic Hamiltonian and the value of the spin. For the 
square-lattice antiferromagnet, close to fe = (0, 0) and (tt, tt) the magnon spectrum 
takes the relativistic forms Bfc = c|fc| and |7r — fe|, c being the spin- wave velocity. 
If we expand n as (l,ei,e2), where the ei are small compared to unity, then the 
equations of motion following from (43) show that there are two modes both of 
which have the dispersion = c|fcj, as expected. If we expand the Lagrangian 
to higher orders in ej, we find that there are interactions between the spin waves 
whose strength is proportional to c/ps, which is of order 1/5. We thus see that all 
the parameters appearing in (43) can be determined by SWT. Compared to the 
standard 1/5* expansion, the liydrodynamic approach is more generic in two points, 
i.e. (i) it is applicable to magnetically disordered phases, and (ii) it may lead to 
non-perturbative results which are beyond the reach of SWT (see, e.g. [51, 52, 53]). 

Ground-state parameters of the 5=1 square-lattice Heisenberg antiferromag- 
net have been studied in great detail using a variety of techniques, including SWT, 
QMC, and series expansions [15]. An early QMC study by Reger and Young [54] 
indicated that the SWT gives a good quantitative description of the ground state. 
Series expansions around the Ising limit performed by Singh [55, 56] found the re- 
sults ps « 0.18J and c « 1.7J, both in good agreement with the first-order SWT 
[6]. Later on, higher-order calculations demonstrated that the second-order correc- 
tions in l/S' to the parameters ps, c and rn are small - even in the extreme quantum 
limit S = ^ - and improve the SWT results. For instance - using both the Dyson- 
Maleev and Holstein-Primakoff formalisms up to second order in l/S - Hamer et 
al. calculated the ground-state energy Ey,/N and the sublattice magnetization m 
[36]. Both formalisms were shown to give identical results closely approximating 
previous series estimates [57]. Different scientific groups have presented consistent 
second-order SWT results for the spin- wave velocity c [58, 59, 60], the uniform 
transverse susceptibility X-l[59, 61] and the spin stiffness constant ps^ [59, 61]. In 
Table 2 we have collected some of these results, demonstrating an excellent agree- 
ment with recent high-precision numerical estimates [37] obtained by using the 
stochastic series expansion QMC method for L x L lattices with L up to 16. 

The accuracy of SWT may be understood in terms of the spin-wave interaction 
V . Indeed, let us consider the 1/5 series for m [36] 

The reported third-order SWT result for this parameter is 0.1750(1)[61]. 
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Table 2. Second-order SWT results for the ground-state energy per site eo = Eq/N 

[36], the on-site magnetization m [36, 59], the spin- wave velocity c [59, 60], the 
uniform transverse susceptibility x± [59, 61], and the spin stiffness constant ps 
[59, 61] of the S = ^ square- lattice Heisenberg antiferromagnet. The SWT results 
are compared to recent stochastic series expansion QMC estimates for Lx L lattices 
with L up to 16 [37]. The series risults for eo, m and x± are taken from [62], and 
those for ps and c - from [61]. The figures in parentheses show the errors in the last 
significant figure, h = ao — J — 1. 



Quantity 


SWT 


QMC 


Series 


-eo 


0.669494(4) 


0.669437(5) 


0.6693(1) 


m 


0.3069(1) 


0.3070(3) 


0.307(1) 


c 


1.66802(3) 


1.673(7) 


1.655(12) 


X± 


0.06291(1) 


0.0625(9) 


0.0659(10) 


Ps 


().18()<)7cS 


0.175(2) 


0.182(5) 



^ = 5- 0.1966019 + + (44) 

For S = ^, the related series in powers of A simply reads m = 0.3033981A" -|- 
0.003464A -|- O(A^), so that the spin-wave interaction V introduces numerically 
small corrections to the leading approximation. The same conclusion is valid for 
the other parameters. 



4.2 Triangular-Lattice Antiferromagnet 

The Heisenberg antiferromagnet on a triangular lattice with nearest-neighbor ex- 
change interactions is a typical example of strongly frustrated spin model.® 

After a long period of intensive studies - see, e.g. [64] and references therein - 

it is now widely accepted that the classical coplanar ground state survives quantum 
fluctuations. This state may bo represented by the ansatz 

= ^ cos(qrM ■ r) + x sin(qM • r) , (45) 

where qm = (^,0) is the wave vector of the magnetic pattern, x _L i: are unit 
coordinate vectors in the spin space, and r runs on the lattice sites. As usual, the 
lattice spacing ao is set to unity. The classical spins lay in the (a;, z) plane, and 
point in three different directions so that the angle ^ is settled between any pair 
of spins in the elementary triangle (sa, Sb, Sc). 

In performing the 1/S expansion about non-coUinear reference states such as 
(45), one faces some novelties which will be discussed in the remainder of this 
Sect. One of them concerns the number of boson fields needed to keep track of 
the whole niagnon spectrum. This is an important practical issue, as higher- order 
spin-wave expansions involving more than two boson fields rule, techni- 

cally intractable. In the general case, this number should be equal to the number 

® For a recent review on frustrated quantum magnets, see [63]. 
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of spins in the magnetic elementary cell, so that for the magnetic structure (45) 
we would need three boson fields. However, in several special cases we can trans- 
form the non-coUinear magnetic structures into a ferromagnetic configuration by 
applying a uniform twist on the coordinate frame. These special systems have the 
property that their magnon spectrum has no gaps at the boundaries of the reduced 
magnetic Brillouin zone connected to the magnetic pattern. The triangular-lattice 
antiferromagnet fulfills this rule, so that we may describe the system by a single 
boson field, as in the ferromagnetic case. In the remainder of this Sect, we shall 
follow this approach [65]. 

To that end, let us rotate the spin coordinate frame about the y axis by the 
angle O^^i — qm ■ [r — r ) for any pair of neighboring spins (s,., s^i ), in accord to 
the reference state (45). In the local reference frame, the Heisenberg Hamiltonian 
acquires the form 

H=Y, [cose„. {sUl'+sUl')+s^^e^r' {sUl'-slsl,)+slsl^ , (46) 

(r ,r ) 

where the sum runs over all pairs of nearest-neighbor sites of the triangular lattice. 

Next, using the Holstein-PrimakofF transformation (4)® and the procedures de- 
scribed in Sect. 2, we find the following boson representation for (46) 



Ti-B = — 



^S'^N + 35 ^ \Aualau + ^ (a+alj, + afea_fe)] + V , (47) 

k 

Ak = 1 + Bfc = —2>uu/2, and Vk = ^[coskx + 2cos(fci:/2) cos(\/3fcy/2)]. Here 

and in the remainder of this Sect., fc takes A'' values from the first Brillouin zone of 
the triangular lattice. 

Up to quartic anharmonic terms, the expansion of the square root in (46) pro- 
duces the following spin- wave interaction V = V3 + V4, where 

V3 = — — y=}{Ki + K2){a\ala3 - a|a2ai) , (48) 



(49) 



' 2 2ViV^^ 

V4 = - X] |-^i2;34«i«2a3a4 + T^'j] (a|4aja4 + ala:ia2ai) 
1-4 

-(1) 



Kfc = i [sin kx-2 sm{kx/2) cos{V3ky /2)] , F^l^^^ = Avis + 4,V2-3 + 1^1 +1^2 + 1^3 + 1^4, 

and _rj23 — — 2(1^1 + !/2 + )• For simplicity, in the last expressions we have omitted 
the Kronecker 5 function, and have used the abbreviations for the wave vectors 
introduced in Sect. 2.2. 

A novelty here is the triple boson interaction V3 = 0{S^^^), which is typical for 
systems exhibiting non-collinear magnetic patterns. We shall see below that such 
kind of interactions complicate the calculation of higher-order 1/S corrections. 



^ The choice of the transformation is a matter of convenience, as the final results 
- at least to second order in 1/S - are independent of the boson representation. 
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Linear Spin Wave Approximation 

In a standard LSWT, we discard V and diagonalize the quadratic part of (47) by 
the Bogohubov transformation ak — Uk{ak — Xka^_k)- The parameters Uk and Xk 
are defined by (16) and (17), but in this case rjh = — 3ffe/(2 + i/k)- The diagonal- 
ization yields the free-quasiparticle Hamiltonian Ho = 3S ixJfea^afc, where the 
dispersion function 

Ek = 3SuJk = 35^(1 -i^fe)(l + 2z/fc) (50) 

gives the magnon energies in a LSWT approximation, to be compared with the 
magnon spectrum resulting from the approach using three boson fields [66] . It is easy 
to check that the dispersion function (50) exhibits three zero modes, as it should be 
since the Hamiltonian symmetry 0(3) is completely broken by the magnetic pattern 
(45). Two of these modes are at the ordering wave vectors k = ±qM, whereas the 
third zero mode at fe = describes soft fluctuations of the total magnetization. 
Expanding about the zero modes, we find the following expressions for the spin- 
wave velocities [67] 

/3\3/2 _ 3^3 „ 

cox=c±g„ = l-j b, Co\\ = Ck=o = . (51) 

Let us now calculate the on-site magnetization m = (s^) = S — (alar). Using 
the Bogoliubov transformation, we find for the density of particles (a|,afe) = — 1/2-|- 
l/(2efe), so that the LSWT result for m reads [66] 

m = S +- - ^ =S- 0.2613 . (52) 

For 5 = i, the LSWT result is m = 0.2387. Since the reported leading 1/5 
correction to m is small and positive''^, there is a clear disagreement with the recent 
QMC estimate m = 0.20(6) [69]. 



Spin Wave Interactions 

Here we consider as an example the calculation of 1/5 corrections to the magnon 
spectrum (50) . There are two different types of corrections related to the spin- wave 
interactions V3 and V4 in (48). Turning to V4, notice that we have already learned 
(Sect. 2.3) that the required correction may be obtained by expressing Vi as a 
sum of normal products of quasipaxticle operators: the diagonal quadratic terms 
give the required 1/5 correction to the spectrum. However, in several cases we 
arc not interested in the quasiparticle representation of V4. Then, it is possible to 
follow another way by decoupling the quartic operator products in V4. Actually, 
this procedure takes into account the so-called one-loop diagrams, and may be 
performed within a formal substitution of the operator products, such as o| 0^0304, 
by the following sum over all the non-zero pair boson correlators 

We are aware of two such calculations reporting, however, somewhat different 
corrections, i.e. 0.0055/5 [68] and 0.00135/5 [65]. 
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a|a2a3a4 i — > ^ ^ [{aJa2)o3a4 + a\a\{a-ia4) — (ala^) {a3a4)'\ ■ (53) 

pair 

As suggested by the quadratic form in (47), there are two types of boson correlators, 
i.e. {a[a2) and (aiag) = (alal), contributing in (53). The constant terms in (53) give 
first-order corrections to the ground state energy, whereas the quadratic operator 
products renormalize the coefficients Ak and Bk in (47). Thus, the interaction V4 
renormalizes the bare dispersion function utk to 

cok = ^JAI - Bl , (54) 
where the new coefficients A]^ and B^ can be expressed in the form* 

^-•^'('+S)+S' ^-«'('+l)+ll- 

An analysis of (54) indicates that the renormalized spectrum still preserves the 

zero mode at fe = 0, but at the same time acquires non-physical gaps at fc = 
±qM. The reason for such kind of behavior of the SWT is connected with the fact 
that we have omitted the l/S corrections resulting from V3. Indeed, the spin-wave 
interaction V3 has the order 0(5^^^), so that a simple power counting indicates 
that IjS corrections to Wfc appear in the second-order of the perturbation theory in 
Va- Wo shall skip the details of this calculation, as it may be performed entirely in 
the framework of the method presented in Sect. 2. Namely, one should express Vij in 
terms of quasiparticle operators, and then apply the general perturbation formula 
(32) for the interaction V3, by using the dressed dispersions (54). As a matter of 
fact, as we are interested in corrections up to l/S, we can use the bare dispersion 
function (50). The final result of this calculation shows that the XjS correction 
resulting from V3 exactly vanishes the gap (produced by V4), so that the structure 
of magnon spectrum (50) - containing three zero modes - is preserved in the leading 
first-order approximation [70]. Based on the renormalized dispersion, the following 
expressions for the spin-wave velocities (51) have been reported [65]: 

0.115\ 0.083\ 

Notice that the l/S corrections diminish the ratio c\\lcx_ from the LSWT result 
1.41 to the value 1.16. These expressions indicate that the leading corrections to the 
magnon spectrum are numerically small even in the case S = \. Good convergence 
has been found also for the IjS series of the magnetic susceptibilities x± and X[| 
[71, 72] which appear as parameters of the magnetic susceptibility tensor [73] 

Xccfi = xx'^c^i + (xii - x^)yc^)yB ■ 

Here y is a unit vector directed perpendicular to the basal (x, z) plane of the planar 
magnetic structure. 

Summarizing, the available SWT results point towards a good convergence of 
the perturbative spin-wave series in the triangular-lattice Heisenberg antiferromag- 
net. This is remarkable, as the spin-wave expansion might have been expected to 
fail for strongly frustrated magnetic systems. 

* For brevity, here we omit the expressions for the constants ai, 02, 61, and 62 [65]. 
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5 Modified Spin Wave Theories 

Here we consider some modifications of tlie standard spin-wave theory allowing 
for an analysis of magnetically disordered phases. These may appear cither as a 
result of quantum fluctuations - a classical example being the spin-S' Heisenberg 
antiferromagnetic chain discussed in Sect. 3.1 - or due to thermal fluctuations, as in 
ID and 2D Heisenberg magnets with short-range isotropic interactions [14]. For the 
antiferromagnetic chain, we have indicated that the failure of SWT arises already 
in the linear spin-wave approximation as a divergency in the boson-occupation 
numbers = (a|ai) = oo implying (sf ) = —oo . Infinite number of spin waves also 
appears at T > 0, when the T = magnetic phases of low-dimensional Heisenberg 
systems do not survive thermal fluctuations. Actually, the occupation rmmbcrs rii 
should not exceed 2S - as dictated by the spin algebra - and the magnetization 
should be zero, as required by the symmetry of the phases. In the remainder of this 
Sect, we discuss modifications of the SWT based on ad hoc constraints imposing 
fixed number of bosons. 

The first generalized spin-wave theory of this kind has been formulated by Taka- 
hashi to study the low-T thermodynamics of ID and 2D Heisenberg ferromagnets 
[74, 75]. Takahashi's idea was to supplement the standard SWT of Heisenberg fer- 
romagnets with the constraint imposing zero ferromagnetic moment at T > 0: 

JV 

M = J2i'''i^)=SN-J2i4ak)=0. (55) 

n=l k 

Depending on the context, in the remainder of this Sect. (A) means the expectation 
value of the operator ^ at T = or T > 0. Quite surprisingly, it was found an 
excellent agreement with the Bethe-ansatz low-temperature expansions of the free 
energy and magnetic susceptibility for the S = ^ Heisenberg ferromagnetic chain. 
Similar extensions of SWT have been suggested for Heisenberg antiferromagnets 
both at T = [76, 77] and at T > [78, 79], by using the same constraint equation 
(55) but for the sublattice magnetization. Below we discuss some applications of 
the modified SWT to low-dimensional Heisenberg antiferromagnets both at T = 
and at finite temperatures. 

5.1 Square-Lattice Antiferromagnet at Finite T 

Using the Dyson-Maleev transformations (9) and (10), the boson Hamiltonian Tig 
of the square-lattice antiferromagnet reads 

'H'B = -YJ^S^+J2\-'^''^^kak+bibk)+Bk{albl+akbk)]+VDM, (56) 

k 

whereas the constraint equation for the total sublattice magnetization takes the 
form 

J2{4<^k+bibk) = SN. (57) 

k 

The wave vector k runs in the small (magnetic) Brillouin zone I kx i ky I < 7r 
containing N/2 points. Ak — JSz, Bk = JSzjk, 7fe = ^{cos kx + cos ky), and a = 4 
is the lattice coordination number. 
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In essence, the constraint equation (57) introduces an effective cut-off for un- 

physical states [80]. To sec this, let us consider the S = ^ system. According to 
(57), the average number of, say, the a magnons is N/4, whereas the total number 
of one-magnon states in the magnetic Brillouin zone is A''/2. Thus, after introduc- 
ing the constraint (57), the effective number of allowed states in the boson Hilbert 
space is 

(N/2)\ 
_(iV/4)!(A74)! 

so that with logarithmic accuracy the correct dimension 2" is restored. 

To implement the constraint equation in the theory, we introduce, as usual, 
a chemical potential fi for the boson fields, i.e. instead of Tig we consider the 
Hamiltonian 

=7^'a-;x^(otafe-h6l6fc), (58) 

k 

where n is fixed by the constraint equation (57). Notice that the introduction of a 
chemical potential simply renormalizes the coefficient Ah — > Ah — /i so that we can 
apply the formalism from Sect. 2 without any changes. 

Using the Bogoliubov transformation (15) with the parameter r)k = JzSfk /{JzS— 
fj,), one finds the following quasiparticle representation oIHb (see, e.g. [17]) 

nB=Eo+nD + VDM, (59) 

where Eo is the ground-state energy calculated up to first-order of the perturbation 
theory in 1/ S: 

E, = -^zJS'(l + l^)\ , = 1-1^0^. (60) 

h 

As we know from Sect. 2.3, the free-quasiparticle Hamiltonian 

Ud = ^Ek{aiak+l3ll3k) (61) 

h 

includes the diagonal quadratic terms resulting from V^jj^, so that the magnon 
energies E^ are calculated up to first-order corrections in 1/5: 

Eu = Jzs{l + ^^ sjl-nl- (62) 

Here the factor r /2S is Oguchi's correction to the magnon spectrum [6]. 

We want to treat the spin- wave interaction up to first order in the 1/5 per- 
turbation theory, so that the Dyson-Maleev interaction Vdm will be discarder. It 
is important to notice that here the off-diagonal quadratic interaction V2 does not 
appear, as dictated by the sublattice interchange symmetry. This means that the 
lowest-order corrections to the sublattice magnetization m have the order 0{S^'^), 
see the series (44), so that the constraint equation (57) calculated in a LSWT 
approximation can be safety used at this level. 

Turning to the magnon spectrum (62), we see that the chemical potential in- 
troduces a spectral gap A so that close to the zone center the excitation spectrum 
acquires the relativistic form 



42^ 

TT N ' 
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E,^^A^+c^k\ c=:^(l + ^), (63) 

where A = 2c{~^/ J zSY^'^ and c is the spin-wave velocity calculated up to first or- 
der in 1/ S. Using the standard expression for free bosons rife = {a\oik) = (Pl-Pk) = 
[exp{—Ek/kBT) — the constraint equation (57) takes the form 

At low T, the main contributions in the last sum come from small wave vectors so 
that, using (63), the gap equation (64) yields 



Zi = I = 2rarcsinh 



1 / 2nps 
oexp f 



A kBTj 



(65) 



Here ps is the T = spin stiffness constant calculated up to first order in 1/S, 
and ^ is the spin correlation length. This result exactly reproduces the saddle-point 
equation in the expansion of the 0(N) nonlinear a model in 2 -|- 1 space- 

time dimensions (see, e.g. [81]). It is well known that (65) describes three different 
regimes, i.e. (i) the renormalized classical, (ii) the quantum critical, and (iii) the 
quantum disordered regimes [53] . 

As an example, we consider the renormalized classical regime defined by the 
condition knT <C Ps. In this case, the last equation yields the following result for 
the correlation length 

(.6) 

This coincides with the one-loop approximation of the 2-1-1 nonlinear a model 
[48]. As is well known, at a two- loop level the T dependence in the pre-exponential 
factor disappears, whereas the exponent argument does not change. 

Finally, let us calculate the leading temperature correction to the internal energy 
U = (Hb)- The expression for U reads 

C/ = i5o + ^^fe (coth£|;-l) . (67) 

k 

Using (63), after some algebra one finds the following result: 

U = Eo+'-^T^. (68) 

Here ("(a;) is the Riemann zeta function. The above temperature correction describes 
the contribution from two zero modes, i.e. k = (0, 0) and k = {it, tt), and reproduces 
the expected universal behavior known from the 2 -|- 1 nonlinear a model and the 
chiral perturbation theory [49, 82]. 



5.2 Applications to Finite-Size Systems 



The modified SWT can also be applied to finite-size systems [76, 77]. This opens 
an opportunity directly to compare SWT results with finite-size numerical data. 
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As is known, the standard SWT is not applicable to finite systems due to diver- 
gences related to the Goldstone zero modes. Actually, the divergency comes from 
the Bogoliubov transformation (15) which is not defined for these modes. 

Turning to the example from Sect. 5.1, notice that in the infinite system the 
chemical potential /j, goes to zero as T ^ 0. At T = the constraint equation takes 
the form 

Here we have selected the contribution from the two zero modes at fe = (0, 0) having 
= ±1. 

According to (69), on a finite lattice the parameter r]o = JzS/{JzS — /i) is 
less than unity, so that the divergences associated with the zero modes disappear. 
The constraint (69) takes into account the fact that in finite systems there are no 
spontaneously broken continuous symmetries. 

To find the staggered magnetization m appearing in the thermodynamic limit 
of the 2D system, we calculate the antiferromagnetic structure factor ^(Tr, n) for 
large N: 

where we have again selected the contribution from the zero modes. 

In the large- TV limit, the last sum transforms into an integral which is oc IniV, 
so that the main contribution comes from the first term in (70). Thus, we find the 
relation 

J™ T^^jTTH- (71) 



2 4 

^ = lim rrr „ 

JV^cx; (1 ■no)N^ 



Equation (69) induces a gap in the magnon spectrum which is defined hy A = 

c-\/2(l — rjfj). Using (71) and the notations from Sect. 5.1, we find the following 
result for the magnon excitation gap in the large- limit 

L = N^^^ is the linear size in a square geometry. The last expression reproduces 
the result for A obtained by other methods [83, 84, 49]. 

Finally, let us return to the Heisenberg antiferromagnetic chain discussed in 
Sec. 3.1, this time using the modified SWT [79]. We have seen that in ID the 
expression for the staggered magnetization (37) contained an infrared divergency 
indicating that the magnetic order is destabilized by quantum fiuctuations. Using 
the concept of the modified theory, we can resolve the problem by replacing (37) 
with the constraint equation 



■ TJq cos^ 



where Kijjo) is the complete elliptic integral of the first kind. 

Since isr(r?o) > 7r/2, the gap equation (73) has a solution for arbitrary S. How- 
ever, the constraint introduces an excitation gap, so that the discussed theory makes 

sense only for integer S. To find the gap, we may use for small (1 — ?7o)^'^^ the asymp- 
totic result Kijio) = In 4(1 — ri'o)~^^^, so that the excitation gap reads 



26 



Nedko B. Ivanov and Diptiman Sen 



Zi ~ c exp {—!tS) . (74) 

Here c is the spin- wave velocity of the antiferromagnetic chain (36). The obtained 
gap has the asymptotic form A ^ 5'cxp(— 7r5), to be compared with Haldanc's re- 
sult A ^ S'^ exp(— ttS*) obtained from the a-model mapping [85, 86]. It is remarkable 
that the simple modified SWT is capble to reproduce the asymptotic expression for 
the Haldane gap. 



6 Concluding Remctrks 

We have surveyed the spin-wave technique and its typical applications to Heisenberg 
magnetic systems in restricted geometries. In most of the cases the SWT results were 
compared with the available numerical estimates. As a result, the systematic large-S 
technique has been found to give very accurate description of the zero-temperature 
parameters and magnon excitation spectra of a number of low-dimensional quan- 
tum spin models, such as the Heisenberg antiferromagnet on square and triangular 
lattices and various quasi-one-dimensional mixed-spin Heisenberg systems exhibit- 
ing ferrimagnetic ground states. Presented analysis of the asymptotic series up to 
second order in the parameter 1/5' implies that in these systems the spin- wave 
interaction introduces numerically small corrections to the principal approxima- 
tion, even in the extreme quantum limit S = ^. Thus, indicated effectiveness of 
the spin-wave technique - as applied to magnetic systems with small spin quantum 
numbers and in restricted geometries - may be attributed to the observed weakness 
of spin-wave interactions. 

The authors thank J. Richter and U. SchoUwock for their collaborations in this 
field, and S. Sachdev, A.W. Sandvik, and Z. Weihong for the permission to use 
their results. This work was supported by the Deutsche Forschungsgemeinschaft. 
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